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We study the inner product of Bethe states in the inhomogeneous periodic XXX spin- 1/2 chain of 
length L, which is given by the Slavnov determinant formula. We show that the inner product of an 
on-shell M-magnon state with a generic M-magnon state is given by the same expression as the inner 
product of a 2M-magnon state with a vacuum descendent. The second inner product is proportional to 
the partition function of the six-vertex model on a rectangular L x 2M grid, with partial domain-wall 
^ boundary conditions. 
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1 Introduction 

The computation of the inner product of two Bethe wave functions is a necessary step on the way of 
obtaining the correlation functions in integrable models. An expansion formula for the scalar product 
of two generic Bethe states in the generalized su(2) model is due to Korepin [1]. The sum formula 
becomes increasingly difficult to tackle when the number of pseudoparticles becomes large. Simpli- 
fications occur when one of the states is on shell, i.e., when its rapidities satisfy the Bethe equations. 
In this case the inner product can be cast in the form of a determinant. A determinant formula for the 
norm of an on shell Bethe state was originally conjectured by Gaudin [2] and later proved by Kore- 
pin [1]. A determinant formula for the scalar product of one on shell and one generic Bethe states was 
found by N. Slavnov [3]. 1 The Slavnov formula was used to obtain some correlation functions in the 
XXX and XXZ spin chains [4], and more recently in the computation of some correlation functions 
in the M = 4 supersymmetric Yang-Mills theory [6,7,9, 10] 

If the two Bethe states are composed from M pseudoparticles with rapidities u = {ui, . . . , um} 
and v = {v\, . . . , vm}, the Slavnov inner product is given, up to a simple factor, by a M x M 
determinant. In [9, 10] one of the authors derived a a closed expression of the inner product for 
the XXX spin chain in the classical limit M — > oo, in which the Bethe roots condense into several 
macroscopic strings. 2 An intriguing feature of this expression is that it has the same form as the 
classical limit of the scalar product of a Bethe state with a vacuum descendent, obtained previously 
in [7], if one chooses the rapidities w of the Bethe state as w = u U v. 

In this paper we prove that the above is true also for finite M. We show that the Slavnov inner 
product (v|u) of two M-magnon states in a periodic inhomogeneous XXX chain with spin 1/2 is 
equal to the inner product of a vacuum descendent with the 2M-magnon state | u U v) , which is in 
turn proportional to the partition function of the six-vertex model with partial domain wall boundary 
conditions, studied recently in [13]. 

2 Inner product of Bethe states in the inhomogeneous XXX 1/2 chain 

The XXX spin chain is characterized by an i?-matrix Ru(u, v) acting in the tensor product V\ <g) Vi 
of two copies of the target space. Up to a scalar factor, the i?-matrix of the XXX spin chain is [14] 



where P\2 is the permutation operator acting in the tensor product V\ <g> V2 of two copies of the target 
space. The inhomogeneous XXX spin chain of length L is characterized by background parameters 
(impurities) z = {zi, . . . , zl} associated with the L sites of the chain. For the spin 1/2 chain, the 
monodromy matrix T a (u) € End(V r a ) represents the product of the i?-matrices along the spin chain, 



The homogeneous XXX spin chain corresponds to the limit z m — > i/2? 

1 Previously, a representation of the inner product as a determinant of double size was found by Kirillov and Smirnov [5]. 
2 The classical limit was studied in [11]; in the condensed matter physics it is known also as Sutherland limit [12]. 
3 Sometimes 0j — Zj — i/2 are referred to as inhomogeneity parameters. With this definitioon the homogeneous limit 
corresponds to 9j — 0. 



Rl2(u, v) = U — V + iP\2 , 



(2.1) 




(2.2) 
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The matrix elements A, B, C, V are operators in the Hilbert space V = V± <8> ■ • • x Vl of the spin 
chain. The commutation relations between the elements of the monodromy matrix are determined by 
the RTT relation 

Rn{u - v)Ti(u)T 2 {v) = T 2 (v)Ti(u) R 12 (u - v), (2.3) 
which follows from the Yang-Baxter equation for R. In components, 

A(v)B(u) = s=2±i B(u)A(v) - ^_B(v)A(u), 

V(v)B{u) = ™ B(u)V(v) + ^B(v)V(u), (2.4) 
[C(u),B(v)} = ^ (A(v)V{u) - A(u)V(v)) , etc. 

As a consequence of (2.3), the families of operators B(u), C(u), as well as the transfer matrices 

T(u) = Tr a [T a {u)} = A{u) + V{u) , (2.5) 

are commuting. 

In the Algebraic Bethe Ansatz, the Hilbert space has the structure of a Fock space generated by the 
action of the creation operators B(u) on the pseudo- vacuum | | L ) = | tt t )> where all spins are 
oriented up. The pseudo-vacuum is an eigenstate for the diagonal elements A and V and is annihilated 
byC: 



(A(u) - A(u)) |t L ) = (£>(«) - D{u)) |t L ) = C(u)\t L ) = 0, (2.6) 



where 



A(u) = Y[ (u - z m + i) , D(u) = [] (u - zm) ■ (2-7) 

m=l m=l 

A creation/annihilation operator with u — > oo is the global su(2) lowering/raising operator, 

«£?Uis-, (2.8, 
L>{u) u A(u) u 

The dual Bethe states are generated by the action of the C-operators on the dual pseudo-vacuum 
(t L | = ( tt t |> which is annihilated by the ^-operators, and the hermitian conjugation can be 
defined as 

C(u) = -B\u). (2.9) 
The space of states is a closure of the linear span of all vectors of the form 

M 

|u) = n^)lt L ), u = {«i,... u M }, M^#u. (2.10) 



The operator B(u) can be viewed as a creation operator of a pseudoparticle (magnon) with momentum 
;h 

generic Bethe states, 4 



p = log "^2 - Such states are called generic, or off shell, Bethe states. The scalar product of two 



4 Here we abuse slightly the established notations, since with the convention B(uf — —C(u) , the state dual to |v) is 
(— l) jv (v|. In our notations the norm of a Bethe state is ||u|| 2 = (— l^^ulu). If one of the states is on shell, then the 
set of its rapidities is invariant upon complex conjugation, and the inner products (— l) JV (v|u) and (v|u) differ only by a 
phase factor. 
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M M 

(u|v) = ^ L \\\C{v k )\\B{u k )\t) i (2.11) 

k=l j=l 

can be computed by applying the relations of the RTT algebra (2.3). For example, the scalar product 
of two one-magnon states is 

( v \u) = (A(v)D(u) - A{u)D{v)) . (2.12) 

A Bethe state is an eigenvector of the transfer matrices [Eq. (2.5)] if the rapidities u = {u± , . . . um} 
satisfies the on-shell condition, which is given by the Bethe equations [14] 

J p-"* + » = _fi " J --^ + » > J = 1,...,M. (2.13) 

Uj U k l llj z m 

k=l J m=l J 



To avoid lengthy formulas, throughout this paper we will use systematically the following notations. 

N 



For any set w = {wj}^ =1 of points in the complex plane, we define the Baxter polynomial 



Qw(«) = Y[(u-Wj), iV = #w, (2.14) 



as well as the rational function 



p±/„.\ 4ef Qu(u ± l) 



E ^-^r- (2 - 15) 

In these notations, the eigenvalue T u (u) of the transfer matrix on the on-shell state |u) is 

T u (u) = Q z (u + i) E-(u) + Q z (u) E+(u). (2.16) 
Another way to write the Bethe equations is as 

e 2i Pu («) = _j ue u (2 .17) 

where the pseudomomentum p u , known also as counting function, is defined modulo tt by 



e 2ipu = 4r — • (2-18) 

As shown by Slavnov [3], when the state |u) is on shell, the inner product with a generic Bethe 
state (v| is a determinant. One can write the Slavnov determinant formula as 

M 

(v|u) = Y[A{ Vj )D{ Uj ) , (2.19) 

j'=i 



det jk £l(uj,v k ) 
deto 



^u,v = ' y f ' , (2.20) 
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where the Slavnov kernel Q(u, v) is defined by 

Q(u, v) = t(u -v)- e 2ip ^ v) t(v - u) , t(u) = — . (2.21) 

u u + i 

For example, Eq. (2.20) gives for the inner product of two one-magnon states (M = 1) 

1 _ E&i. 

(v\u) = A{v)D{u) y U)V = i . (2.22) 

This expression indeed matches with the restriction of the general expression (2.12) when the rapidity 
u is taken on shell, D(u)/A(u) = 1. 

3 Alternative expression for the inner product 

3.1 Operator factorization formulas 

The Slavnov determinant (2.20) can be given a very convenient operator expression [9, 10], whose 
derivation we review below. We represent the Slavnov kernel Q(u, v) as the result of the action of two 
difference operators on the Cauchy kernel u ^ v+i , 

Q(u, v) = (1 - e 2 ip»(v) e id / dv ) ( e - id / du - l) 1 (3.1) 

and write the Slavnov determinant as the result of the action of N pairs of difference operators to the 
Cauchy determinant , 

& _ Ht'gvl 1 K e ) lineup 1 ) ut3L jfc uj-vk+i 

=>u,v — i • yj- 1 ) 

ael jk Uj -V k +i 

Here, and in the following, the formulae contain products of difference operators and the ordering of 
the difference operators should be respected. The factors within each of the two blocks in the above 
formula commute, but the factors belonging to different blocks do not. 
Now we apply the Cauchy identity 



det 1 = X\ 1<k ( u 3 ~ u k) Uj<k( v k-Vj) = A u A-y 

jk Uj -v k + i Uj!k=i( u j -Vk + i) ~ n u,v 

After repeated application of the obvious identities 



e -id/du 1 = 1 e - l9 /du ( u G u) 

' = E-{ Vj )-^-e id ' d ^ (v€v) (3.4) 



TT -U K-J/ TJ 



and taking into account the expression (2.18) for the pseudomomentum, we write Eq. (3.2) in a 
factorized operator form, 

y u , v = (-!)#» ±- [J (l - §^ e»/»>) A v ■ ± J] (l - Et(u) e-»l°») A u . (3.5) 
Here, we have to be careful in that the operator exp(id/dv) acts on all factors on the right of it. 
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3.2 The functional #/ u [f] 

The two blocks of factors in the above operator expression have a similar form and suggest introducing 
the following quantity. For any set of points u = {uj}jL 1 in the complex plane and for any complex 
function f(z), we define the functional 

^ = ^-n( i -/(«) e±<o/ou ) A - 0.6) 

Substituting A u = detjk(uj~ l ), one can write this functional as a ratio of determinants [9, 10] 

<[f] = — rTA ( 3 - 7 ) 



det 



jk 



but for our purposes the operator representation (3.6) is more convenient. 

The functional stf^ [/] can be expanded as a sum of monomials associated with the partitions of 
the set u into two disjoint subsets, 

<m= e <«» n [-/(«% <u-= n e u>( u ")- (3.8) 

u'Uu"=u «"Gu" «"6u" 

Under this form, the functional £?^[f] appeared previously in Ref. [7]. If the function f(u) depends 
implicitly on u, we define the functional si^ [f] so that it is given by the same expansion (3.8). 
The operators srf+ [/] and stf~ [f] are related by the functional identities 

<-[/] = <[v/] nt-zwi' <w =<"[v/] n [-/(»)]' (3 - 9) 

ti£u nSu 

<-[/] = <[-g/], <+[/]=<-[- |^/] • (3-10) 



The first pair of identities was proved in [9, 10]. Here we give the proof of the second pair. 
ProofofEq.(3.10): We transform the coefficients srf~, u „ of the expansion (3.8) of s>/~[f] as 

<-,„»= n *»-(«-)= n PS £ >"' = ( - 1)#u " n PS w>. < 3i » 

where we used the property = E^,E^„, as well as the obvious identity 



n PSr ( - 1)#u '' am 

Summing over all partitions, we obtain the expansion of the functional £/ + [g], with g = —E~/E+ f. 
□ 
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3.3 A symmetric expression for the inner product 

We can associate with the functional s/^ 1 [/] a difference operator s^^ [/] acting on the functions on 
the set u, by replacing f(u) — > /(«) e ±t9 / 9u for all u € u. The operator =e^j t [/] is well defined if 
the function f(u) does not depend implicitly on the variables u. The c-functional =e/ u [/] is the result 
of the action of the operator functional [/] on the constant function 1, 

< ± [/]=< ± [/]-l- 

Then the expression [Eq. (3.5)] for the inner product can be laid down in terms of the functionals srf~ 
and as 

^u,v = (-l) #u K + [E+/E+] ■ s/-[E+] . (3.13) 

Eq. (3.13) is equivalent to the operator factorization formula derived in [9, 10]. 

Below we give an alternative expression for the inner product, symmetric with respect to the sets 
of rapidities u and v. Define the functional 

^ u ,v = (-l)^*t+[E+/E+]-sS+[E+/E+] =y v , u . (3.14) 

If the rapidities u in Eq. (3.13) are assumed to be on shell, we can apply the first identity (3.10) to 
transform the second factor on the rhs as 

,3 = 0) Bf\ *£>> <[-e^§] (2 i 7 » <[§]. (3.15) 

Therefore we can identify 

y u v = y u v for u on shell. (3.16) 

The functional ^ u ,v can be simplified further. It is given, up to a sign, by the functional £/ w [1 / E+\ , 
with w = u U v: 

^u,v = (-l) #u *Cj v [1/E+]. (3.17) 
Proof: Using the definition (3.6) and decomposing 

A uUv = A u A v H (u-v), (3.18) 

u£u,t)£v 

we write s/^j v \\/E^\ as a product of two operators, one depending explicitly only on the set u, and 
the other depending on the set v: 

-£ n (i - §g ^) ^ ■ £ n (> - E M e ' 8/ i A - 

= J v + [e+/e+] • *C[e+/e+] (3 - 19) 
= (-i) #u ^u,v 

As an example, apply this formula for the inner product [Eq. (2.12)] of two one-magnon states: 

= .«£,,) [1/-EJ] =1 - E+(u)/E+(u) - E+(v)/Et(v) + l/£ z + (t.)E z » 
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If the rapidity u is on shell, then is = 1 and the first term disappears, while the second gives the 
inner product. 

A representation of the rhs of (3.17) as a Fock expectation value of chiral fermions is given in the 
Appendix. 

4 The Slavnov determinant as a pDWPF 

The inner product of two M-magnon states can be thought of as as a partition function of the six- vertex 
model on a 2M x L rectangular grid, with particular boundary conditions. The i?-matrix R(u, v) [Eq. 
(2.1)] can be graphically represented as two intersecting segments, a horizontal one carrying a rapidity 
u and a vertical one carrying rapidity v, 



K(u - v) = " Z 



j 


k 




m 



= {u- v)5 jk 5l m + i Sjl5km . (4. 1) 



The two halves of each segment are labeled by su{2) indices, which are represented in a standard 
way by arrows. Then the six non-vanishing elements of the i?-matrix correspond to the vertices of the 
six- vertex model in the rational limit [16], depicted in Fig. 1, with Boltzmann weights 

a(u, v) = u — v + i, b(u,v)=u — v, c(u,v)=i. (4.2) 

The operators B(u) and C(u) are graphically represented in Fig. 2. The inner product [Eq. (2.11)] can 
be identified as the partition function of the six-vertex model defined on a 2M x L rectangular grid, 
shown in Fig. 3. The partition function is a sum over all possible ways to associate arrows with the 
internal links, so that at each site of the lattice the number of the incoming arrows equals the number 
of the outgoing arrows. The boundary conditions on the arrow configurations are the following: on 
the two vertical boundaries, the lower half of the arrows point outwards, while the upper half point 
inwards. On the two horizontal boundaries all arrows point upwards. For example, the scalar product 
of two one-magnon states, Eq. (2.12), is represented as 



(v\u) = 



k=i 





L J 















































(4.3) 



Z L-1 Z L 



The identity (3.16) means that the Slavnov inner product has a second interpretation as a partition 
function of the six-vertex model. Before formulating the correspondence, we give a brief recollection 
of the (partial) domain wall partition functions. 

The domain wall partition function, DWPF, is the partition function of the six-vertex model on a 
square grid with domain wall boundary conditions, DWBC [1, 17]. The DWBC restrict the arrows 
on the edges forming the two horizontal boundaries to point outwards, while the arrow on the edges 
forming the two vertical boundaries point inwards, as shown in Fig. 4. The weights depend on 
the rapidity variables w = {wi, . . . , wjy} and z = {zi, . . . , zn}, associated respectively with the 
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Figure 1 : Graphical representation of the 
non- vanishing elements of the i?-matrix rep- 
resented by the vertices of the six-vertex 
model. 



Figure 2: Graphical representation of the 
elements of the monodromy matrix 



Zj z 2 



Ui - 

U 2 - 

"m- 
Vl - 
V 2 - 



z, z 



1 L 2 



w 2 



W N - 



Figure 3: Graphical representation of the Figure 4: The domain wall boundary con- 

inner product (v|u) as a six-vertex partition ditions (DWBC). 

function on a rectangular grid. 



z, Z 2 



w 2 



z, z 



1 L 2 



u 2 

U M 
Vl 

y 2 



Figure 5: The partial domain wall bound- 
ary conditions (pDWBC). The boundary ar- 
rows on the top point outwards, those on the 
vertical boundaries point inwards, and those 
on the bottom segment are free. 



Figure 6: The pDWBC relevant for the 
Slavnov inner product. 



horizontal and with the vertical lines. As shown by Izergin [18, 19], the DWPF can be expressed, up 
to a factor which can be eliminated by a renormalization of the six- vertex weights, as a determinant 5 



_ tot^-H) 1 1 

det,i. — it n + i 

The Gaudin-Izergin determinant (4.4) is a particular case M = L = N, u = w, v = z, of the 
Slavnov determinant (2.20). Indeed, if v = z, the second term in the Slavnov matrix SI vanishes, 
Sl(u, v) = t(u — v). Then the first factor in (3.5) is equal to 1, and we find [10] 



'W,Z 



■l) N */-[E+] = (-l) N s/+[E-], #w = #z = N. (4.5) 



Eq. (4.5) can be derived also directly from the definition (4.4), proceeding in the same way as we 
did in the case of the Slavnov determinant. In this derivation one does not need to assume that the 
rapidities w are on shell. 

The quantities named in [13] partial domain wall partition functions, or pDWPF, are defined on 
a rectangular grid n x N, where 1 < n < N. The boundary conditions are the same as the DWBC 
except for the bottom segment of the boundary. On the N x N square grid, the positions of the bottom 
arrows are fixed, once we have imposed the DWBC on the left, the right and the top segments of the 
boundary. When n < N, this is no more the case. The partial DWBC, depicted in Fig. 5, do not 
impose any restriction to the positions of the arrows on the bottom segment, and the bottom arrow 
configurations are summed over, just as those for the bulk segments. We will denote the pDWPF 
again by 2f Wj2 :, with 1 < #w < #z. When this is needed by the context, we will denote by w ra a set 
of rapidities w with cardinality #w = n. 

The pDWPF can be obtained from the DWPF by sending the rapidities w n+ i, . . . , wn, associated 
with the bottom N — n rows, sequentially to infinity. Let u n be the subset of the first n rapidities in 
ujv = {u±, . . . , un}- Then the result of sending the remaining N — n rapidities to infinity is [10, 13] 

= u^o T "' uJSoo T^™ = ( " 1)n(iV - ^ ^6) 

Applying the second identity (3.9) to the rhs of (3.17), we can relate the functional <^ u ,v to the 
partial domain wall partition function i? u uv,z defined on the rectangular 2A^ x L grid, as the one 
shown in Fig. 6, where N = #u and L = #z, 

^u,v = (-l) #u §^ , AA w , z d ^ f (#z-#w)!n^ + H. (4-7) 

Another way to write the proportionality factor is as 

A/w, z = (#z-#w)! [J T)t\' < 4 - 8 ) 

The proof Eq. (4.7) is surprisingly simple. First we transform the representation (4.6), using the 
first of the functional identities (3.9), 

iT uUv , z = (L — 2M) ! <- Uv [E+] = A/- uUv , z <+ Uv [1/E+] . (4.9) 

This concludes the proof of (4.7). 



5 For the first time the ratio of determinants (4.4) appeared in the works of M. Gaudin [2, 20] as the scalar product of 
two Bethe wave functions for a Bose gas with point-like interaction on an infinite line. 
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5 Concluding remarks 



In this paper, we derived an expression for the inner product of an M-magnon Bethe eigenstate and 
an M-magnon generic state in the inhomogeneous periodic XXX chain of length #z = L, which is 
completely symmetric in the union of the two sets of rapidity parameters: 

^ u ,v = (-l) #u < + Uv [l/£+]. (5.1) 

The functional S^ u ,v represents a determinant 2M x 2M and is given essentially by the partition 
function with domain boundary conditions on a L x L square grid, with L — 2M of rapidities sent to 
infinity. The functional ^ UiV coincides with the original Slavnov product (3.13) if the rapidities u are 
on shell, and for periodic boundary conditions (no twist in the Bethe equations). In general, ^ U)V and 
y u , v are two distinct functionals. 

An immediate application of the representation (5.1) in the study of the semiclassical limit of 
the three -point function of long trace operators in the su(2) sector N = 4 SYM, formulated in 
Refs. [6,7]. Using the determinant representation given in [8], a closed expression for the structure 
constant for three non-protected operators was obtained in Refs. [9, 10] as a generalization of the result 
for one-protected and two non-protected operators operators found in Ref. [7]. The classical limit of 
the functionals y U:V and ^ u ,v is the same, but if one is interested in the subleading terms, second 
functional is much more convenient to deal with. On the other hand, the expression (5.1) with generic 
inhomogeneity parameters can be used to reproduce the higher orders in the weak coupling expansion 
of the structure constant, as it has been argued in [30, 31]. 

The alternative representation of the inner product (2.11), found in this paper, has a natural in- 
terpretation on terms of the Fock space for the Algebraic Bethe Ansatz. The functional ^ u ,v is 
proportional to the inner product 

M M 

^ U ,v ~ (i L i(s-) L - 2M n^)n^)it L ) 

3=1 3=1 

M M 

~ (t L \ HC(vj) Y[C(uj) (S+) L - 2M U L > . (5.2) 

3=1 3=1 

The second functional has the same structure as the result of a particle-hole transformation on the 
ket vector on the rhs of the original inner product [Eq. (2.11)]. In case of a non-zero twist, the dual 
rapidities are different than the original rapidities [29]. In the case we are considering, M of the dual 
rapidities coincide with the original ones, while the the rest L — 2M of them go to infinity. This can be 
justified by the following simple argument. 6 For zero twist, the global su{2) symmetry is not broken 
and the Bethe states is a direct sum of states with given spin, which are eigenspaces of the transfer 
matrix. The states corresponding to the same solution of the Bethe equations, belong to the same 
su(2) multiplet [32]. Therefore the Bethe eigenstates Y[ ueu B{u)\ | L ) and ~\\ u( z u C{u)\ ], L ) must be 
related by the action of a global raising operator, which must be (§ + ) i_2M , since the first state has 
S z = \L- M, while the second state has S z = M - \L. Thus we have 

M M 

Y[C(uj) (8 + ) L - 2M \i L ) ~ Y[B( Uj ) \t L ), (5.3) 

3=1 3=1 



6 We thank D. Serban for suggesting to us this argument. 
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where the proportionality sign means that the two states are equal up a c-function of the rapidities. In 
this sense, our main result is the computation of the factor of proportionality. 



Note added. After the preparation of the manuscript, we learned about the publication [33], which is 
along the same lines as our work. 
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A Fermionic representation of the functionals srf^- [/] 

The functionals £/ u {E+] and £/+[l/E+] have natural fermionic representations. Introduce a 
chiral Neveu-Schwarz fermion living in the rapidity complex plane and having mode expansion 

lp(u) = ^2 1prU~ r ~^, 1p(u) = ^2 VV" r ~5. (A.l) 

The fermion modes are assumed to satisfy the anticommutation relations 

[VV, 4>r'] + = [VV, 1pr']+ = , [A, i>A+ = 5 r,r> , (A.2) 

and the left/right vacuum states are defined by 

(0\ip- r = (0|Vv = and ip r |0> = ^_ r |0) = 0, for r > 0. (A.3) 

The operator ip r creates a particle (or annihilates a hole) with mode number r and the operator ip r 
annihilates a particle (or creates a hole) with mode number r. The particles carry charge 1, while the 
holes carry charge —1. The charge zero vacuum states (A.3) are obtained by filling the Dirac sea up 
to level zero. The left vacuum states with integer charge ±N are constructed as 

f (o|y>! ...Yv_i if iv > o, 

(N\ = { _ 5 (A.4) 

Any correlation function of the operators (A.l) is a determinant of two-point correlators 

(0|^(«)^(«)|0) = ^. (A.5) 
The following formulae is easily established. Let #u = n and #z = N. Then 

= ^-niny =l[ ^)-^ + o] hslm) 

(N-n\UU^)Uti^k)\0) 
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^_ rz?+1 {N-n\Y[ n j=1 ty{ Uj )-^{u j + i)] Uti^kM 

£?,, \EJ\ = jrf — = . (A.7) 

(N-n\UU^ j + i)Uk=i^)\0) 

Eq. (A.6) gives a convenient representation of the pDWPF and the Slavnov scalar product in 
terms of free chiral fermions. Other fermionic representations have been proposed in [21-23] for the 
Slavnov product and in [24-26] for the Gaudin-Izergin determinant, see also the review paper [27]. 
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